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- Abstract 

o ; 

. A Charge Transport (CT) mechanism has been proposed in several papers {e.g., Yavin et al. PNAS 102 

' 3546 (2005)) to explain the localization of Base Excision Repair (BER) enzymes to lesions on DNA. The CT 

. mechanism relies on redox reactions of iron-sulfur cofactors that modify the enzyme's binding affinity. These 

' redox reactions are mediated by the DNA strand and involve the exchange of electrons between BER enzymes 
along DNA. We propose a mathematical model that incorporates enzyme binding/unbinding, electron transport, 
and enzyme diffusion along DNA. Analysis of our model within a range of parameter values suggests that the redox 

. reactions can increase desorption of BER enzymes not already bound to their targets, allowing the enzymes to be 

qh' recycled, thus accelerating the overall search process. This acceleration mechanism is most effective when enzyme 
I— copy numbers and enzyme diffusivity along the DNA are small. Under such conditions, we find that CT BER 

^ ' enzymes find their targets more quickly than simple "passive" enzymes that simply attach to the DNA without 

^ . desorbing. 

' Key words: Repair Enzymes; Charge Transport; Target Search; Lesion; DNA 
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(N ■ 

m ; 1 Introduction 
in ■ 

' The genomes of all living things can be damaged by ionizing radiation and oxidative stress. These factors can cause 
0> I mismatches in the DNA strand resulting in localized lesions. The role of Base Excision Repair (BER) enzymes 
' is to locate and remove these lesions. If the lesions arc allowed to persist, they can give rise to mutations and 
^ , ultimately diseases such as cancer, 
r ■ The localization of BER enzymes to lesions is physically related to the binding of transcription factors to 

rS I promoter regions that regulate gene expression. In 1970, experiments by Riggs et al. P, 0] showed that the 
' association rate of the Lad repressor to its operator is about 100 times faster than the maximum rate predicted 
by Debye-Smoluchowski theory. This theory assumes that Lad is transported to its target on DNA via 3D 
diffusion. To explain the experimental observations, the theory was modified to account for facilitated diffusion 
0, 0, 01- In this process, the Lad repressor can spend part of its time attached to the DNA and perform a ID 
random walk before detaching and diffusing in 3D again (see Fig. [T]). Provided the protein spends about half its 
time on the DNA and half its time in solution and the diffusivities in ID and 3D are comparable, the predicted 
search time can be reduced by as much as 100 fold [sj. However, these conditions are very restrictive as the 
protein can spend up to 99.99% of its time associated to the DNA 0] and the diffusion constant along DNA (in 
ID) is in general much smaller than the one in the cytoplasm (in 3D) Therefore, many modifications of the 
basic facilitated diffusion theory have been proposed, including intersegmental transfers [9|], the effect of DNA 
conformation [lo| , directed sliding [ll| . and finite protein concentration |l2| . 

A series of recent papers 1^ 14l. Iisj have revealed a special kind of long-ranged interaction for certain BER 
enzymes based on charge transport (CT) along DNA. MutY, a type of DNA glycosylase, contains a [4Fc-4S]^+ 
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cluster which is very sensitive to changes in its environment. Specifically, its redox potential is modified depending 
on whether it is in a polar environment (when the enzyme is in solution) or in a more hydrophobic one (when the 
enzyme is attached to DNA). In solution, the [4Fe-4S]^+ cluster is resistant to oxidation. However, when attached 
to DNA, the cluster is more easily oxidized through the reaction [4Fe-4S]^+ ^ j4Fe-4S]'^"'" + e~. Furthermore, the 
3+ form has a binding affinity about 10,000 times greater than the 2+ form [l5|. 

A model for the "scanning" of BER enzymes along DNA, aided by CT, was proposed in 1^, 14, 15|, and is 
depicted in Fig. [2l When a BER enzyme adsorbs to DNA, it oxidizes and releases an electron along the strand 
(see Fig. Eta))- Distal enzymes, already adsorbed onto the DNA can absorb these electrons, become reduced and 
desorb. Hence, binding and unbinding of enzymes are associated with oxidation and reduction of their iron-sulfur 
clusters. CT along DNA can be disrupted by the presence of defects that affect electron transport. For example, 
guanine radicals ("oxoGs"), formed under oxidative stress, can absorb electrons: see Fig. [l^b). By acting as 
sites of reduction, they promote the adsorption of BER enzymes 13, [l^. Once the radical has absorbed an 



electron, it converts to a normal guanine base and no longer participates in CT. However, "permanent" defects, 
or lesions, can also exist on DNA which can absorb more than a single electron (see Fig. [SJc)). For example, 
oxoGs can erroneously pair with adenine bases when the DNA replicates. Such lesions may continuously absorb 
electrons with a certain probability, or otherwise refiect them. In contrast to the oxoG-cytosine case, the removal 
of oxoG-adenine lesions require MutY to be present at the damaged site. 

In this paper, we develop a model of CT- mediated BER enzyme kinetics that includes enzyme diffusion along 
DNA, a binding rate that depends on electron dynamics, and the effects of finite enzyme copy number. Our 
key finding is that the proposed charge transfer mechanism employed by BER enzymes accelerates their search 
for targets along DNA in real finite enzyme copy number systems. In the next section, we derive the governing 
equations of enzyme kinetics. These equations are rendered non-dimensional and key parameters are defined and 
estimated. In Section 3, we numerically solve our model equations under various conditions and estimate the time 
for the binding of an enzyme to a localized lesion. We conclude with a discussion of our results in Section 4. 



2 Mathematical model 

2.1 Derivation of kinetic equations 

Consider the diffusion and adsorption-desorption kinetics of repair enzymes in a bacterium such as E. coli: see 
Fig. [TJ The chromosome in bacteria is circular but tightly coiled up into a nucleoid that has an effective volume 
of about 8 X 10^ nm^. If a repair enzyme is associated with the DNA strand, it can diffuse freely along the DNA 
to find lesions. These associated enzymes can spontaneously desorb from the strand, but they can also become 
oxidized, leading to tighter binding to the DNA. If later on, the enzyme is reduced, its association with the DNA 
weakens and it can quickly dissociate from the DNA. Localized lesions prevent the passage of electrons (released 
along the DNA by oxidation of associated repair enzymes) by either reflecting or absorbing them. 

We write mass-action equations for the reactions occuring in Fig. [3l coupled to equations that determine the 
electron dynamics. We assume that the enzyme density in the bulk, Rb{t) (where t is time), is well mixed and 
has no spatial dependence. The density of DNA-adsorbed BER enzymes in the reduced and oxidized state are 
denoted by Ra{x,t) and Q{x,t) respectively, where < a; < L is the coordinate along the DNA and lesions are 
located at a; = and x — L. The density of guanine radicals is g{x,t) and the density of rightward and leftward 
electrons is N+{x,t) and N-{x,t). Note that Rb{t) has units of inverse volume, while Ra{x,t), Q{x,t), N±{x,t), 
and g{x, t) carry units of inverse length. The governing equations corresponding to the processes depicted in Figs. 
[2] and [3] are 
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= D^^~v{N^ + N^)Q + mR^, (1) 



dt dx^ ' ^ " \L 



L 



-J^ = -Konitfc + / Radx, (3) 

aAMM)^„aiV^ . + (4) 

afcii = (6) 

These equations must be solved subject to the boundary conditions 



(7) 



(8) 



N+{0,t) = riV_(0,i), N^{L,t) ^ rN+{L,t), 
Q(0, t) ^ Q{L, t) = 0, i?,(0, t) = Ra{L, t) = 0, 

and initial conditions 

Q(a;,0) = 0, i?a(x,0)=0, i?b(0) = no/17, 
7V+(x,0)=0, iV_(x,0) = 0, g(x,0) = go/i. 

In Eqs. [TJini D+ is the diffusivity of adsorbed MutY3+ along the DNA, £»_ is the difFusivity of adsorbed MutY^^, 
V is the speed of electrons along DNA, m is the electron release (oxidation) rate of adsorbed MutY^^', fcoff is the 
intrinsic desorption rate of MutY^"'' , fcon is the intrinsic adsorption rate of MutY^"*" to the DNA from solution, O 
is the cell volume, L is the arclength of the DNA, and / is the electron flip rate (see below). In Eqs. [7l r is the 
electron reflectivity of lesions, which we describe in more detail later. In Eqs. [3 no is the copy number of MutY, 
and go is the initial number of guanine radicals on the DNA. The definitions of all constants arc summarized in 
Table [H 

We now give a brief justification of equations [l][6] and conditions [7] and [H The form of the first three equations 
can be understood from Fig. [Hb) which summarizes the reactions among the three species Ra, and Q. Eq. 
[T] describes the time rate of change of adsorbed MutY'^^ due to oxidation of adsorbed MutY^^ (+mRa) and 
reduction by incoming electrons {—v{N^ + N-)Q). The first term on the right hand side represents diffusion along 
the DNA. Eq. [2]describes the evolution of adsorbed MutY^+ in terms of the reduction of MutY^+ {+v{N++N^)Q), 
spontaneous desorption into solution (— A:offi?a), adsorption of aqueous MutY^^ {koa{Q/L)Rb) and oxidation into 
MutY3+ i-mRa). Since MutY^^ binds to DNA less strongly than MutY"^'*', it is possible that D+ is appreciably 
smaller than D_. Eq. [3] is an equation for the concentration of MutY^^ in solution which can decrease by enzymes 
binding to the DNA {—ko^Ri,) and increase by enzymes unbinding from the DNA (represented by the integral 
term). Because we assume enzymes in the bulk solution are well mixed, any increases in bulk concentration are 
due to an integrated DNA-adsorbed density which does not distinguish between enzymes that are released from 
different positions along the DNA, but only sees the total number of enzymes that desorb. 

Eqs. |4] and [5] describe the electron dynamics. In our model, right and left-moving electrons (sec Fig. [31^ a)) 
propagate along the DNA with speed v; this process is represented by the two convectivc terms on each of the left 
hand sides. Also, electrons are lost when they are absorbed by MutY^"*" or by guanine radicals, and produced when 
released by adsorbed MutY^^. These processes are represented by the third and fourth terms on the right hand 
side of Eqs. m and m respectively. Finally, leftward and rightward electrons can inter-convert by scattering off 



inhomogeneities and thermally induced conformational changes in the DNA [25|, |26| . This process is represented 
by the flrst and second terms on the right hand side. The flip rate / characterizes how frequently a traveling 
electron changes direction. If / is large, the electron move diffusively, but if / is small, it moves in a more ballistic 
manner. Finally, Eq. [6] represents the evolution of the guanine radical population. OxoGs arc annihilated when 



^The mean free path of an electron is estimated to be A ~ 1 — 10 base pairs and the flip rate approximated as v/X. 
^Estimated using the time taken for the restriction endonuclease BsoBI to unbind from DNA [23|, toff =150s and taking fcoff = 1/ 
This value of tos may not be an accurate value for the unbinding time for MutY. 

^Assumes that about 1 in 40,000 guanine bases are oxoGs [2^] and the length of E. coli DNA is L = 5 x lO'' bp. 
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Symbol 


Definition 


Typical value 


Reference 


D± 


Diffusivity of adsorbed enzymes 


5 X 10*^ bpVs 




V 


Electron velocity 


10^° bp/s 


m 


f 


Electron flip rate 


109 - lOio s-i 


1 


m 


Electron release rate 


~ 10^ s^i 


[19] 


n 


BacterimxL volume 


3.7 X lO^nm^ 


[20] 


L 


Length of DNA 


5 X 10^ bp 




h 

"'Oil 


MutY^"^ attachment rate 


2000 s^i 


[16] 


Ks 


MutY^"'" detachment rate 


7 X 10~3 s"i 


2 


no 


Copy number of MutY in E. coli 


20-30 


[21, 22] 


r 


Electron reflectivity of lesions 


- 1 




90 


Number of oxoGs on E. coli DNA 


30 


3 



Table 1: Key constants for used for repair enzyme model Eqs. [T]-[U]and conditions 



Parameter 


Deflnition 


Calculated Value 


rj 




< 10-^° 


U 


vj (konL) 


5 


a 


m/{m + kos) 


~ 1 


F 


f /kon 


5 X 10^ - 5 X 10^ 



Table 2: Dimensionless parameters in Eqs. [TOlfni 



Accelerated Search of DNA Repair Enzymes 



5 



they absorb electrons as represented by the —v[N^ + NJ)g term. Radicals might also be spontaneously generated 
and modeled by a source term on the right hand side of Eq. [6l In this paper, we neglect spontaneous oxoGs 
generation. 

Equations [T][5] use a mean field approximation that neglects stochastic fluctuations in enzyme, electron and 
guanine number. The effect of noise in the system could be included through the use of a chemical master equation 
[27[ ; however, generalizing the equation to account for spatial variations along the DNA is beyond the scope of 
this paper [28|. Nonetheless, we expect our results for lesion targeting by enzymes will be qualitatively accurate. 

Lesions at a; = and x = L (see Fig. El^a)) define the domain of solution for Eqs. [T][6] which are subject to the 
boundary conditions [7l This domain can represent a circular DNA with circumference L containing a single lesion. 
In the first equation of [71 leftward traveling electrons are converted to rightward traveling ones by the lesion that 
reflects leftward electrons with probability r. If r = 0, leftward electrons are absorbed by the lesion. On the other 
hand, if r = 1, the lesion is fully reflective and the rightward and leftward electron densities arc equal. Similar 
considerations apply to the lesion a,i x = L. Since we will eventually use our mean- field mass action equations to 
estimate the mean time for a repair enzyme to find a lesion, we assume that the lesions are perfectly "absorbing" 
and set Q ~ Ra = at the lesion positions. Our simulations are performed on a domain with go oxoG radicals 
and a bulk solution that contains hq enzymes (see Eqs.[8|); hence the adsorbed oxoG density is go/L and the bulk 
concentration is uq/Q. 

2.2 Model reduction and non-dimensionalization 

Before non-dimensionalizing Eqs.[T][6]wc can make one important simplification. On the right hand side of Eq. [21 
the sizes of the second, third, fourth and fifth terms arc approximately v/L'^, kog/L, k^^J L and m/ L (in units of 
bp^^ s^^) respectively. Guided by Table[Tl we assume the term mRa dominates. Since the oxidation rate is large, 
adsorbed MutY^"'' quickly oxidizes into the 3+ form upon adsorption onto DNA. For times t ^ 1/m and rates 
fcoff + v/L, koa Eq. [2] gives Ra{x, t) for all < x < L and we can neglect spatial gradients in Ra as well 
as dRa/dt. Therefore, we approximate Eq. [2lwith 

Ra{x, t) « i— (v{N+ + iV_)Q + fcon f ■ (9) 

Upon substitution of Eq. [9]into Eqs. [H [4] and [3 we eliminate the equations for Ra and find a reduction analogous 
to one commonly used in deriving the steady-state limit of Michaclis-Menten kinetics [2^ . 

We now non-dimensionalize our equations by measuring time in units of k~^l , length in units of L, concentration 
of adsorbed species in units of 1/L and concentration of bulk species in units of 1/17. Our final set of reduced 
and nondimensionalized equations that describe the transport and kinetics of MutY repair enzymes, right and 
left-moving electrons, and guanine radicals is 

^^^"^'^^ - -U{l-o){N++N^)Q + 7^i^ + nR,, (10) 



dt V /V + . -j^ , 

dRbit) 



dt 



1 

U{1-<t) I {N+ + N^)Qdx - (jRb, (11) 



dN,ix,t) , ^dN^ ^ FiN.-N,)^gUN,+ ^-^-(l-^^)uN,Q+^^UN.Q, (12) 



dt dx ' "^^ ^ 2 V 2/ 2 

dN.[x.t) ^mjx^ ^ -FiN^-N,)-gUN^ + ^-^ + ^-UN,Q-(l-^-)uN^Q, (13) 



at dx ' ^' ■' 2 2 \ 2 

^ .u{N+ + N^)g, (14) 
where we have defined the dimensionlcss c^uantities 

"Tin i^nu J-' "Tin 
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and 

m 

TO + fcoff ' 

which can be estimated using Table [21 As we discuss later, the parameter a represents the effective binding rate 
in terms of the competition between the electron release rate m and the dcsorption rate of DNA-bound MutY^+ 
fcoff, and lies between and 1. The dimcnsionless boundary and initial conditions are 

iV+(0,t) = riV_(0,t), 7V-(l,t) = riV+(l,i), Q(0, t) = 0(1, t) = 0, (17) 

and 

Q(.T,0) = 0, Rb{0) = no, , . 

N+{x,Q) = N_{x,Q) = 0, 5(x,0)=5o. ^ ^ 

Our model can approximate the case of "infinite" enzyme copy number when the transport of bulk enzymes is 
diffusion-limited. Although most of the enzymes cannot immediately adsorb onto the DNA as they are too far 
away, we assume that a certain number, Rb, are in the vicinity of the nucleoid, say within a volume CI' (see Fig. 
[T]), and are able to directly engage in adsorption. However, instead of being depleted over time, Rb is continuously 
replenished by far enzymes that diffuse into fi' C to keep Rb fixed. Therefore, to obtain the infinite copy number 
limit, we hold Rb constant in Eqs. [TOl [121 and [T31 and Eq. [Til no longer applies. To summarize, we model the 
infinite copy number case by holding Rb constant. In the finite copy number case, Rb{t) is allowed to vary in time 
through Eq. [TTl Finally, note that equations describing a simple diffusing enzyme that does not undergo CT can 
be recovered from Eqs. [TOKMI bv setting U = 0. In this case, the equations for Q{x,t) and Rb{t) decouple from 
the rest. 



2.3 Repair enzyme binding affinity cr 

In Eqs. [TmiTil the rate of creation of reduced, adsorbed enzyme Ra from reduced bulk enzyme Rb is exactly Rb 
since we measure time in units of 1/fcon- However, the overall rate of the compound reaction Rb ^ Ra — > Q is 
aRb- Consider a MutY"^^ that is adsorbed onto the DNA. If it absorbs an incoming electron, it can either desorb 
into the bulk or it can release an electron back along the DNA and remain oxidized. The parameter a in Eq. 1161 is 
the probability of electron release. When fcoff ^ m, a MutY'^^ that absorbs an electron will preferentially desorb 
{Ra Rb), but when fcoff <C to a MutY"^^ will simply release the electron it just absorbed to stay adsorbed onto 
the DNA {Ra ^ Q)- These limiting behaviors are realized by taking cr ^ and a 1 respectively. 

If (T 0, a bulk reduced enzyme that adsorbs onto the DNA quickly desorbs back into the bulk, while if cr 1, 
MutY^^ on the DNA prefers to oxidize and stay adsorbed rather than go into solution. Once it is oxidized, any 
further electrons that are absorbed will be re-emitted in a random direction. Hence the electron changes direction 
with probability 1/2 whenever it encounters an adsorbed MutY"^^: when cr = 1, the terms with prefactors (1 — cr/2) 
and cr/2 in Eqs. [T^ and [Ol add to the F{N^ — terms to yield an effective flip rate of F + UQ/2. The seeding 
of oxidized enzymes on the DNA increases the effective electron flipping rate because these enzymes can absorb 
electrons and immediately release them back along the DNA in the direction they came from or in the direction 
they were going. 

We end this section with the comment that the model for the CT redox process in Fig. [2lis not exactly equivalent 
to the reaction scheme in Fig. ^h). In Fig[21 a bulk MutY^+ {Rb) adsorbs onto a DNA and immediately oxidizes, 
releasing an electron along the DNA. DNA-bound MutY^+ {Q) remains adsorbed until it absorbs an incoming 
electron, whereupon it reduces and immediately desorbs into the bulk. For this model to hold, the reaction kinetics 
in Fig. [HJb) must be non-Markovian. Specifically, consider the intermediate quantity Ra in Fig. [Sjb). An Ra 
enzyme oxidizes to a Q enzyme {Ra Q) only if it 'remembered' that it was originally created via a, Rb Ra 
reaction. Likewise, an Ra enzyme desorbs {Ra — > Rb) only if it 'remembered' that it was originally created through 
& Q Ra reaction. 



3 Results and Discussion 

We now compute and analyze solutions to Eqs. 11011141 for the infinite and finite copy number cases. The equations 
are solved numerically using second order finite differences on a non-uniform grid that clusters grid points near 
the boundaries and a trapezoidal rule to approximate the integrals. MATLAB's stiff solver odelBs was used to 
integrate the equations in time. In the infinite case Rb is held at the value tiq and in the finite case, Rb{t) is 
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included in the dynamics with initial condition i?b(0) = uq. Furthermore in each case we consider the dynamics 
associated with CT enzymes where U > 0, and the dynamics associated with "passive" , non-CT enzymes where 
U = 0. Setting U ^ decouples the equations for electron and guanine radical dynamics (Eqs. [12lfT4|) from the 
equation for Q{x,t), the density of DNA-bound enzymes (Eq. fTU]) . 

We shall explore the behavior of Eqs. [TUlfHl and the associated search times defined below, with respect to: 

• <T, the effective binding affinity. Generally we have < a < 1. From the values of m and fcoff in Table [1] we 
have a « 10~^. This value of <t renders the desorption term —U{1 — a){N+ + N-)Q in Eq. [TOl insignificant, 
making the effect of CT negligible. Therefore, a necessary requirement for an effective CT mechanism is that 
a < 1. In our simulations for the MutY system, we take a = 0.9 bearing in mind that the value of fcoff in 
Table [U is for BsoBI and not MutY. 

• 77, the diffusivity of MutY^+ along DNA. The value in Table [2] of rj = 10^^° is based on the diffusive sliding 
of a human glycosylase, hOggl, which has a diffusivity of about 5 x 10® bp^/s [l^]. However, this value may 
not necessarily be an accurate value for MutY. Therefore we will explore a range of diffusivities 77 near 10"^". 

• go, the initial guanine radical density: There are about 30 oxoGs at any given time on E. coli DNA, but this 
number depends on environmental conditions. Hence we explore a range of values centered around go = 30. 

• r, the lesion reflectivity: the interaction between electrons and lesions depends on unknown molecular factors 
at the lesion and in the bulk cytoplasm. Hence, we explore a full range of r values between and 1. 

• F, the electron flip rate: the precise dynamics of electrons on DNA is a very complicated process; our 
estimate for F in Table [5] makes many simplifications and may not be accurate. We will explore a range of 
F values centered around 10^. 

Figure |4] shows snapshots of adsorbed enzyme, guanine and electron density profiles for a finite enzyme copy 
number (no = 30) system. The profiles are shown near the lesion at a; = at times t — 2 and t = 5. The 
electron density is generally smaller at the lesions and larger in the middle of the domain, resulting in a larger 
enzyme desorption rate away from lesions (the desorption rate in Eq. [TO] is proportional to the total electron 
density N+ +7V_). Thus, the CT enzyme density is smaller than that for passive enzymes away from lesions. The 
enhanced desorption of CT enzymes from the interior continuously replenishes the number of enzymes in solution 
so that Rb{t) decreases less rapidly than for passive enzymes. For intermediate times, the net deposition rate is 
larger for CT enzymes, the enzyme density near the lesions is also larger (Fig. (Ua)) and grows in time (Fig. [IJb)). 
For long times, the density vanishes everywhere: this is the trivial steady state solution to Eqs. [TOlfHl 

Figure [5] shows the DNA-bound enzyme density a.t t ~ 40. In (a), there is a sharp spike in the enzyme density 
near the lesion at a; = 0, but otherwise the enzyme density is relatively small. Note that all densities are symmetric 
about X = 1/2. In Eq. [TOl CT enzymes desorb with a rate proportional to the total electron density iV+ + N-. As 
seen in (a), this density is smallest at the lesions. Therefore the enzyme density near x — and x = I grows more 
quickly compared to the interior density. The inset shows a rapid variation in Q of about 600 within a boundary 
layer of width ~ 10"'^. Using a non-uniform grid that clusters the mesh points near the boundaries, we are able 
to resolve these boundary layers to calculate the flux of enzymes through the lesions. In Fig. [^b), CT {U = 1) 
and passive (U = 0) enzyme densities are compared when the copy number is flnite. The CT-enzyme density has 
sharp maxima near the lesions, while the passive enzyme density docs not. Compared to the inflnitc copy number 
case, the size of the maxima is smaller since the number of enzymes in the bulk (and hence the deposition rate) 
decreases with time. Because of the maxima, the flux of enzymes into the lesion. 



Jit) = r; 



dQix,t) _ dRaix,t) 



dx dx 



(19) 



x=0 



is greater compared to the passive case. Figure [SJe) compares the current for CT ([/ = 1) and passive enzymes 
{U = 0) when the copy number is infinite. The passive enzyme current is always greater than the CT enzyme 
current because for a constant deposition rate, any desorption reduces the number of enzymes on the DNA and the 
flux of enzymes into the lesion. Therefore, for infinite copy number systems, search by passive enzymes will always 
be faster than CT enzymes. In contrast, when the bulk contains a finite number of enzymes. Fig. EJd) shows that 
the CT current is always greater than the passive current. This is due to free electrons on the DNA that determine 
the local desorption rate. In the CT mechanism, enzymes are knocked off the DNA by incoming electrons and on 
average, desorb from lesion-free portions of the DNA and re-adsorb at positions closer to the lesion. The result is 
that for intermediate times [t > 10), the current experiences a second growth phase, a behavior that is never seen 
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for the passive case. Ultimately, we have J{t) ^ as t ^ oo for both passive and CT enzymes, but CT ensures 
that this behavior occurs at a much later time. 

Next, we consider the typical time for the first enzyme to reach a lesion. Since the enzyme density is symmetric 
about X = 1/2, the total flux can be found by using twice the enzyme flux to one lesion defined in Eq. 1191 The 
typical search time is then approximated by integrating 2J[t) until one enzyme has diffused into the lesion: 



From solving the full set of equations [T][B] numerically, we find that the gradients in Ra{x, t) at the lesions are neg- 
ligible compared to those of Q{x,t), verifying the validity of eliminating Ra and using J{t) ?» f]{dQ{x,t)/dx)x=o 
as the total enzyme current. In the mean field limit, an alternative definition of the search time is Tg ~ 



tive differences from computed using Eq. 1201 

Figure[6l^a) shows that the search times are extremely sensitive to the initial number of oxoGs go. In particular, 
there is a rapid increase in as go increases past the enzyme copy number uq — 30. The CT mechanism relies on 
the presence of free electrons that cause enzymes to desorb from lesion-free portions of the strand and re-adsorb 
near lesion sites, while oxoGs suppress CT by absorbing free electrons. When go > riQ, all enzymes from the bulk 
adsorb onto the DNA and any released electrons are absorbed by nearby guanine radicals. Instead of participating 
in CT-mediated redistribution and localization, the enzymes cannot desorb and must rely on slow diffusive sliding 
along the DNA strand to find their targets. When go < '^Oi 8-1 least one enzyme is always in solution and is 
transported through the cytoplasm. Since 3D transport is assumed to be fast, the search time is correspondingly 
small. Also in this plot, increases as r increases but the search time is much more sensitive to go- the search 
time changes by about 20% for go ~ and by about 0.05% for go ~ 50 over the whole range of r. In our model, 
the search time is not greatly affected by whether lesions refiect or absorb electrons; what is important is that the 
lesions prevent their passage along the DNA. 

Figure [6ljb) shows that for the range of ?/ values explored, there is a value < cr* < 1 for which the search 
time Ts is minimum. To understand why there is an optimal cr*, consider the CT mechanism's dependence on 
the binding affinity a. If cr = 1, enzymes strongly bind onto the DNA. Even when they absorb electrons, they 
will re-emit them to stay adsorbed on the strand. Hence, there is no desorption, no fast transport through the 
cytoplasm and acceleration of the search. On the other hand, if ct = 0, enzymes do not stay on the DNA long 
enough to even slide into lesions and the search is correspondingly slow. Our results show that the search is 
optimal when the enzyme is weakly associated with the DNA i.e. < cr ^ 1 and the electron release rate is small 
compared to the intrinsic desorption rate. From the data in Tables [T] and [21 it seems that real cells do not operate 
near this optimal regime. In Fig. [Hlb), we computed most of the search times using unrealistically large values of 
r] to clearly show the minimum with respect to cr. For smaller rj values we have cr* — > 0+, but the dependence of 
Ts on cr does not change qualitatively. For larger rj values, enzymes do not rely on CT to localize to lesions and 
can find their targets quickly using diffusive sliding. In this case, the search is most efficient if as many enzymes 
as possible adsorb on the DNA; this situation is realized by taking cr = 1 and r, monotonically increases as cr gets 
smaller. 

Figure [Zta) shows how the search time varies as a function of ID enzyme diffusivity along the DNA. Notice 
that the search time for CT enzymes {U = 1) is much smaller than that for passive enzymes (U = 0). Indeed, 
can be reduced by several orders of magnitude when the effects of CT are included. If fewer oxoGs are initially 
present, the search occurs more quickly. Consistent with Fig. Efa), the search time is extremely sensitive to the 
initial number of guanine radicals on the DNA. For passive enzymes, Ts scales as 0{r]~^). For CT enzymes, the 
0{r]~^) behavior switches to Ts = 0{ri~^^^) for sufficiently large ?/ with the crossover dependent on go. 

For finite copy number Fig. [3b) again shows that the search time decreases if CT is included, but this time 
for different ffip rates. For the large values of F used in Fig. [71 one can show that the effective diffusion coefficient 
of the electron density scales as 1/F [l^. Therefore, as F increases the electron density dissipates more slowly 
through the partially absorbing lesions. A greater density of free electrons implies more enzyme desorption, more 
transport through the cytoplasm and faster search times. In the F oo limit, we expect the enzymes to self- 
desorb independently of the oxoG density. This can be seen from Eqs. 1121 and [T3l where the dominant terms on 
the right hand side are ±F{N^ — N-^) and aRb/2. In principle, as _F ^ oo, one can approximate N± in terms of 
Rb{t) and substitute into —U(l — (t){N^ + N^)Q in Eq. [10] to further reduce the system to only two equations 
for Q{x,t) and Rb{t). 




(20) 



t exp — Jg* J{t')dt' dt. We have computed Tg using this mean field approximation and find negligible qualita- 
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Although both plots in Fig. [7] arc for the finite copy number case, we also performed analogous simulations 
for the infinite copy number limit. We found that including the effects of CT by taking U ~ 1 always led to an 
increase in the search time compared to the passive case: for fixed rj and increasing U always increased Tg 
regardless of the value of go . 



4 Conclusions 

Our key finding is that charge-transport (CT) dynamics mediated by redox reactions can significantly reduce 
search times of repair enzymes in real cells where the copy number is finite and the diffusivity along the DNA is 
small. In theoretical systems where the copy number is infinite, CT actually slows down the search. The speed-up 
in finite systems arises because of a spatially dependent desorption rate. Specifically, the desorption is greater 
along intact portions of the DNA but smaller near lesions. As a result, CT-induced enzyme-enzyme interactions 
"recycle" enzymes so that they desorb from lesion-free parts of the DNA and reattach closer to lesion sites. Our 
proposed mechanism is illustrated in Fig. 8. A related mechanism has been implicated in mRNA translation 
where ribosomes are recycled, enhancing protein production rates, (soj 

If we re-dimensionalizing the search times by using an estimated value of fcon — 2000 s~^, we find that passive 
enzymes with diffusivity rj ~ 10"^" have long search times of about 15 minutes (see Fig. [7]), comparable to the 
the life cycle time of E. coli. With the CT mechanism and go = 28 initial oxoGs, the search time drops to a few 
seconds. For smaller values of 77, the difference in search times between passive and CT enzymes becomes greater. 
Using go = 20, we calculate to be about 30 hours and 2 seconds respectively for passive and CT enzymes 
that have diffusivity 77 ^ 10^^^. Therefore for realistic enzyme diffusivities, we think that CT is an indispensable 
mechanism that allows enzymes such as MutY to locate lesions on the DNA in a reasonable amount of time. 

When the initial number of oxoGs exceeded the enzyme copy number, we found a large increase in the search 
time. In this case, search takes place mostly through slow diffusive sliding along the DNA. However, when the 
copy number (number of potential electron emitters) exceeds the number of electron absorbers, we find that the 
search time decreased drastically, with the search taking place mainly through the transport of enzymes through 
the cytoplasm. Therefore, we predict that the spontaneous generation of electron absorbing defects (such as oxoG) 
would significantly slow down the search and conversely, the presence of other redox-active proteins (such as the 
transcription factor SoxR [3l| ) would speed up the search. Although such proteins may not be directly involved in 
lesion search, they may be upregulated when the cell is oxidatively stressed, increasing the population of electron 
emitters in the system. The iron-sulfur cluster responsible for CT in MutY is also found in other repair enzymes 
like EndoIII [s^l- Hence EndoIII could also participate in CT, emit electrons to promote the desorption of MutY 
and speed up the search. 

Recall that classical facilitated diffusion theory [^, 0, 0| predicts a large reduction in the search time of proteins 
providing equal amounts of time arc spent in ID and 3D. However, most proteins are strongly associated with 
DNA so that the speed up is not achieved in practice. The passive enzyme system considered in this study can be 
thought of as a suboptimal search by facilitated diffusion: with U = 0, a MutY that oxidizes and binds to the DNA 
cannot desorb back into the cytoplasm and the protein spends much more time diffusing in ID. However, when 
U > 0, bound oxidized MutY can be knocked off the strand by electrons. CT therefore provides a mechanism for 
MutY to spend more time in 3D than it otherwise would. In other words, CT-aided MutY could be one system 
where the conditions required for speed up are actually satisfied. In addition, when MutY binds near lesions, it 
may diffusively slide along the DNA into its target: the target size is effectively increased with the DNA acting 
like an antenna [lo| . This antenna effect is enhanced by enzymes preferentially oxidizing and adsorbing onto parts 
of the DNA that are near lesion sites. 

Extensions to our model may include spatial gradients in the bulk enzyme concentration, more careful treatment 
of electron dynamics, and adding fiuctuations in copy number. Nonetheless, our simple deterministic model 



describes mechanisms and yields results qualitatively consistent with findings in jl3l. Il4| 
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Figure Captions 

Figure 1: (a) Target search on prokaryotic DNA, which is tightly coiled up into a nucleoid. Proteins in the bulk can 
diffuse to the DNA through the cytoplasm to locate their targets, (b) Searching proteins (hexagons) locate targets 
(diamonds) by sliding along DNA, punctuated by attachment and detachment. Q represents the cell volume while 
Q' represents all points in the vicinity of the nucleoid. Enzymes within Q' can engage in direct adsorption onto 
the DNA. 



Figure 2: Charge Transport (CT) mechanism proposed in 1^, 14, 15|. (a) A repair enzyme (in solution) is in the 
2+ state and adsorbs onto the DNA. Its iron-sulfur cluster oxidizes in the process, releasing an electron along the 
DNA. A repair enzyme (already adsorbed on the DNA) is in the 3+ state and accepts an incoming electron. Its 
iron-sulfur cluster reduces and the enzyme desorbs. (b) Guanine radicals ("oxoGs") can absorb free electrons on 
the DNA. These radicals are annihilated upon absorbing an electron, (c) Lesions can partially reflect and absorb 
electrons. 



Figure 3: (a) Summary of the CT model, described by Eqs. [T][H1 Bulk enzymes, with density Rb{t), can attach to 
the DNA and oxidize to release rightward and leftward electrons with densities N+{x, t) and N-{x, t) respectively. 
Guanine radicals with density g(x, t) act as electron absorbers. Upon adsorption, oxidized enzymes with density 
Q{x,t) are formed with Ra{x,t) as a transient, intermediate quantity. Fixed lesions are located at a; = and 
X = L. (b) Redox reaction diagram for the MutY repair enzyme. MutY^+ in solution is represented by Rb{t), 
MutY^"*" adsorbed onto DNA is represented by Ra and MutY'^+ adsorbed onto DNA is represented by Q. 



Figure 4: Density profiles for enzyme, guanine, and electrons on DNA in a finite enzyme copy number system 
(i?f,(0) = no = 30) at time (a) t = 2 and (b) t = 5. Dashed lines show density profiles of passive enzymes in which 
the CT mechanism is absent. Parameters used were rj = 10^^", a = 0.9, F = 10^, and go ~ 28. 



Figure 5: Enzyme profiles and currents for infinite ((a) and (c)) and finite ((b) and (d)) copy numbers. In each 
figure, the profile or current is plotted for passive (dashed) enzymes where U = and CT (solid) enzymes where 
U = 1. Insets show the large gradients in enzyme density within a thin boundary layer near the lesions. Parameters 
used were cr = 0.9, ?/ = lO^i", .go = 28, r = 0.5, F = 10^ and no = 30. 



Figure 6: Search time for Charge- Transport Enzymes, for copy number no = 30, electron flip rate F = 10^ 
and electron speed U ~ 1. The actual search time in seconds can be recovered by dividing by fcon, whose value is 
estimated in Table [TJ (a) Search time as a function of initial guanine density go and lesion electron reflectivity r. 
Parameters used were a = 0.9 and 77 = 10"^°. (b) Search time Ts as a function of enzyme binding affinity a and 
enzyme diffusivity along DNA, rj. Parameters used were go = 30 and r = 0.5. 



Figure 7: (a) Search time of passive enzymes {U = 0) compared with CT enzymes (U > 0) as a function of 
enzyme diffusivity i] for various initial guanine densities go- Parameters used were a ~ 0.9, r ~ 0.5, F ~ 10^ 
and no = 30. (b) Search time of passive enzymes compared with CT enzymes for different electron flip rates F. 
Parameters used were a = 0.9, r = 0.5, 77 = 10^^° and no = 30. For both plots, the actual search time in seconds 
can be recovered by dividing by k^n, whose value is estimated in Tabled] 
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Figure 8: Recyling of enzymes via the Charge- Transport mechanism. In a finite copy number system, the mecha- 
nism increases the enzyme desorption rate for intact portions of DNA but decreases it near lesions. Therefore on 
average, enzymes are "recycled" to lesion sites by 3D transport through the cytoplasm. In many cell systems, this 
method of finding lesions is faster than a ID search by diffusive sliding. 
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